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Abstract 

We study a relativistic diffusion equation on the Riemannian phase 
space defined by Franchi and Le Jan. We discuss stochastic Ito (Langevin) 
differential equations as a perturbation by noise of the geodesic equation. 
We show that the expectation value of the angular momentum and the 
energy grow exponentially fast. We discuss drifts leading to an equilib- 
rium. As an example we consider a particle in de Sitter universe. It 
is shown that the relativistic diffusion of momentum in de Sitter space is 
the same as the relativistic diffusion on the Minkowski mass-shell with the 
temperature proportional to the de Sitter radius. We study a diffusion 
process with a drift corresponding to the Jiittner or quantum equilibrium 
distributions. We show that such a diffusion has a bounded expectation 
value of angular momentum and energy. The energy and the angular 
momentum tend exponentially fast to their equilibrium values. 

1 Introduction 

In this paper we discuss relativistic dynamics of a particle in general relativity 
which moves in a medium (gas) of other particles. We assume that the interac- 
tion with a gas consists of frequent elastic collisions which in a Markov limit (no 
memory) can be approximated by a diffusion process. In the framework of the 
special theory of relativity the diffusion problem has been formulated and solved 
a long time ago by Schay [1] and Dudley [2] under the assumption that the dif- 
fusion is evolving in the proper time on the phase space and preserves the mass 
. The diffusion problem on a general pseudoriemannian manifold has been 
formulated recently in a geometric framework of diffusions on fiber bundles by 
Franchi and Le Jan [3] . The theory has been further developed and applied in [4] 
[5] [6] [7] (for a general coordinate dependent discussion of relativistic diffusions 
on the pseudoriemannian manifolds see [8]). 
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We follow the method of our earlier papers [9] [10] [11] to define the relativistic 
Brownian motion on the submanifold of fixed of the Riemannian phase space. 
With our choice of coordinates the generator of the process and the stochastic 
equations coincide with those of Franchi and Le Jan [3] . The diffusion process 
describes a random perturbation of the geodesic motion. We show that in a 
spherically symmetric static background the expectation value of the angular 
momentum and the energy grow exponentially fast. The model is rather un- 
physical without a friction. Following our earlier paper [9] we discuss the friction 
terms leading to the Jiittner and quantum distributions (diffiisions with Jiittner 
equilibrium distribution are also discussed in [12] [7]). We show that such a fric- 
tion ensures finite expectation values of energy and angular momentum. 

The diffusion process describes an evolution of the particle in a gas around 
a star (or a black hole ) before achieving an equilibrium with the surrounding 
medium (this may be the gas of Hawking radiation [13] from the black hole). In 
ref.[10] we studied photon diffusion in an electron gas. Here, we discuss a soluble 
example of a particle diffusion in de Sitter universe. We obtain a surprising 
result that the diffusion of momentum is described by the same equation as 
the relativistic diffusion on the Minkowski mass-shell with a friction leading to 
the Jiittner distribution. The corresponding temperature is proportional to the 
radius of de Sitter space. The relativistic diffusion (without friction) in the 
Schwarzschild metric has been studied earlier in [3] and in the Godel universe 
in [4]. For reviews on the relativistic diffusion see [14] [15]. 

The paper is organized as follows. In sec. 2 we discuss the geometry of the 
mass-shell. In sec. 3 we construct the generator of the diffusion on the mass-shell 
as the Laplace-Beltrami operator on the level surface in the cotangent bundle. 
In sec. 4 we restrict ourselves to isotropic metrics. We discuss the stochastic Ito 
equations in these coordinates. In sec. 5 we obtain the transport equations in the 
laboratory time corresponding to the diffusion equations in the proper time. We 
discuss friction forces which lead to equilibrium measures determined by basic 
principles of the statistical physics. In sec. 6 a time cvohition of the energy and 
the angular momentum is discussed. In sec. 7 we take the limit of zero mass. 
The diffusion in de Sitter space is discussed in sec. 8. In two appendices we give 
a derivation of some results applying methods of the stochastic calculus. 



2 Riemannian geometry of the mass-shell 

Wc are interested in random pertiirbations of the dynamics of relativistic par- 
ticles of mass m moving in a gravitational field. The dynamics on a pseudorie- 
manian manifold M. (the geodesic motion) is described by the equations [16] 
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where /k = 0, 1, 2, 3 and 

dsl = g^u{x)dx^dx'" (3) 

is the pseudoriemannian metric onM. The four- momentum p(t) of a relativistic 
particle defines the mass by the relation 

= g'^^'p^.P. = rr?c^. (4) 

If eq.(4) is satisfied then from eq.(l) it follows that r has the meaning of the 

proper time (which is invariant under coordinate transformations). Eq.(4) de- 
fines a submanifold in the cotangent bundle. We take the pseudoriemanian 
metric on the cotangent bundle M. x TM* as the product metric on }A and 
TM* 

ds^ = ds\ + ds\ = gf_t^{x)dx^'dx" + gf"" {x)dpf,dp„. (5) 

The metric on the submanifold (4) is inherited from the product metric on the 
cotangent bundle. It can be obtained by expressing by the spatial components 

of the momenta (here g^"^ is the inverse matrix to the one defining the metric on 
TM; we choose the convention that the spatial terms in ds"^ have an opposite 
sign) 

dsl ^ g^''{x)dpi_,dpu = g'^'^dpodpo + 2g'^''dpodpk- g^'^dpjdpk = -^^''dpjdpk (6) 

{j, k = 1,2, 3) .We assume that 3°° > and g^° = (static metrics can be chosen 
in this form). Then, we obtain 

y fe = gjk _ oj-^gjrgknp^p^ (7^ 

We have 

D = det(7J'=) = det(c/^'=)w-2 (8) 

and 

Ijk = 9jk + m.-'^c-'^pjPk. (9) 

The measure invariant under all coordinate transformations of the subman- 
ifold (4) of TM* reads [17] 

PQ^dxdp. (10) 



3 Diffusion on the cotangent bundle 

We define the Laplace-Beltrami operator on the cotangent bundle as usual by 
means of the metric. The metric (5) is a product of the metrics whereas the 
Laplace-Beltrami operator A is a sum A^^ + A^ . The part of the Laplace- 
Beltrami containing derivatives over the momenta reads 
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In order to define the Laplace-Beltrami operator on the submanifold (4) of the 
cotangent bundle we have to restrict the operator (11) to this submanifold. A 
way to do it is to treat the condition (4) as a definition of the level surface 
in the cotangent bundle. In such a case the Laplace-Beltrami operator on the 
level surface is the same as the one defined by the Riemannian geometry of the 
Riemannian mass-shell discussed in sec. 2. 

AH = i?-^^7,fei?^^, (12) 
dpj dpk 

where D is defined in eq.(8) and jjk in eq.(7). 
Explicitly, 

Ah = {gjk + m-^c-^PjPk)^4-^ + Sm-^c-^Pk^. (13) 



The diffusion equation is 
with the generator 



dr(l)r = a<^r (14) 



(15) 



The operator (15) generates a diffusive dynamics which is a perturbation of the 
geodesic motion (l)-(2) by a diffusion in the momentum space. A coordinate 
free frame bundle definition of the diffusion (14) has been obtained earlier in 
[3](Lemma 3.1). 

We still transform eq.(15) in order to express the diffusion on the Riemannian 
mass-shell as a perturbation of a diffusion on the Minkowski mass-shell (we apply 
this transformation in sec. 8). For this purpose let us introduce the tetrads f" 
by the equation 

= (16) 

(a = 1, 2, 3). We define the inverse fl 

fifj=€ (17) 

We change coordinates on the cotangent bundle 

P, = ffPa (18) 

while = x'^ . The coordinate vector fields take the form 

— = P— (19) 
dp, ^'^dp'J ^''^ 

d d d 
a^ = a^+'-U^, (20) 
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where 

^la = ^fl (21) 

The vector field of eq.(20) can be considered as a horizontal lift of the coordinate 
vector field from M to M x TM*. 



4 Relativistic diffusion in isotropic coordinates 

We shall consider highly symmetric manifolds. Such manifolds admit special 
coordinate systems. We restrict ourselves to the isotropic spherically symmetric 
metric[16] from now on. Then, 

dsl = A-2(|x|)c2df2 _ B-2(|x|)dx2. (22) 

In these coordinates 

Ah = {B-%k + m-^c-^p^pk)^^^ + 3m-^c-^pk^^. (23) 
The diffusion equation in the isotropic coordinates reads 

1 3 In A , ,2 d<l> , K.^m^c'^ a j y^^^) 

where in the generator (15) we expressed po as 

PQ = A-^{B'^Y>'^+m^c^)^=A-^uj. (25) 

We can express the solution of the diffusion equation (24) as an expectation 
value E[..] over the sample paths of a diffusion process {xriPr) starting from 
(a;,p) [18] 

(j)r{x,v)=E[(l){Xr,Vr)]. (26) 

The stochastic process can be defined as a solution of a stochastic equation. In 
order to write down the stochastic equations we must calculate the square root 
of the matrix (7^7), i.e., to find a tetrad e 

e?e?=7,7. (27) 

We obtain 

e« = B-i(5« + — (w - mc)^-^PiPa). (28) 
■' ■' mc 

Then, Ito stochastic differential equations [18] corresponding to the represen- 
tation (26) of the solution of the diffusion equation (24) read (see Appendix 
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A;stochastic equations for the relativistic diffusion liave been obtained also in 
the frame bundle formalism of ref. [3], Theorem 3.2) 



—Auj (29) 
dr m 



rfa;° _ 1 
dr rr 

dx 1 

= B^p 30 

dr m, 

dpj = ^Bd^Bp^dT - ^u^d^ In AdT + ^Pjdr .3^^ 
+mcKB~^dbj + kB~^{lo — mc)p~'^pjpdh. 

Here, b(s) denotes the Brownian motion defined as the Gaussian process with 
values in and the covariance 

E[ba{s)bc{T)] = 6acmin{s,T) 

5 Invariant measure and kinetic equations 

We define an evolution of the measure da = dxdp^ by the equality (see [9]) 

{(l)r)cr = J dxdp^(j)r = J dxdp<^r4>- (32) 

Then, 

dr<^r - Q*^r (33) 

where Q* is the adjoint of Q in L'^{dxdp). We say that the measure a is invariant 
if 

{Cl>r)a = {4>)a (34) 

is independent of r. In such a case 

= 0. (35) 

Eq.(35) is a transport equation in the laboratory time x'^ . A normalizable 
invariant measure needs a dissipation. We add such a dissipation term Kjdr to 
the stochastic equation (31). This is equivalent to adding the friction 



K = K,— (36) 



to the diffusion equation (24). 

In isotropic coordinates eq.(35) reads 

m ^ dx^* ^ ^ m dx^ ^ ^3 ^ 2m dx^ dpj ^ ^ ~^ 2m dx^ dpj ^0 ^ 

(37) 
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In the search of the equihbrium distribution we demand that the dynamic 
part and the diffusion part of the equation (37) cancel separately. So, for the 
dynamic part 

-^A'Po^^E + ^IjB'pj^e - ^^^P'^E + = 0. 

(38) 

Then, from the diffusion part 

Kk = K^m^c^^-^'[^^{B-^d,k+m-^c-^PjPk) - |(mc)-V)^£. (39) 

Any function of deterministic constants of motion is a solution of the dynamic 
equation (38) (static metric with a spherical symmetry has the angular momen- 
tum and the energy as constants of motion). Let us consider the solution as a 
function of po (25) in the form ( Jiittner [19] and quantum statistical distribu- 
tions are of this form) 

= Auj-'^ exp{f{pcA-^uj)). (40) 

From eq.(39) we obtain 

Kk = ^PkHcA-^ujfipcA-^oj). (41) 

(diffusions with Jiittner equilibrium distribution have been discussed also in 
[12] [7]). Explicitly, with the friction (39) the diffusion equation is 

' ^ m ax^' ' rn ^ J ax^ " 2m ax^ ^ cJpj 2m ox^ opj 

+ ^^B-%,+m-^c-'pjPk)^^^ + K'{l + '^ujcPA-^f{l3cajA-^))p^^J. 

(42) 

Then, the stochastic equation for the process Pr (26) solving the diffusion equa- 
tion (42) reads (see Appendix A) 

dpi = ^Bd^Bp^dr - ^to'^dnnAdr 

m 3 ^ m 3 

+ ^ojcl3A-'^fdT + ^PjdT + mcKB-^dbj + kB-'^{oj - mc)p-^pjpdh. 

(43) 

The equilibrium distributions (40) (Jiittner or quantum) are not invariant under 
Lorentz transformations even in the flat (Minkowski) case. The notion of the 
equilibrium depends on the Lorentz frame where the particle equilibrates. We 
make this frame dependence explicit in [20]. 

6 Evolution of the energy and the angular mo- 
mentum 

In the isotropic coordinates we define the angular momentum 

L = X X p. (44) 
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If ^(|x|) and -B(|x|) depend only on |x| then the metric is sphericahy symmet- 
ric. In such a case when k = then the angular momentum is a constant of 
motion. First, we consider the proper time evolution of the angular momentum 
in a stochastic model without a friction. From eq.(31) we obtain a stochastic 
equation for the time evolution of the angular momentum when k ^ (see 
Appendix B) 

dLr = —^Lrdr + kB~^{lo - mc)p~'^Lripdh) + KmcB~^x x dh. (45) 

In the spherical coordinates (r, 6, </>) eq.(45) can be expressed in the form 

3^2 

dp<i, = —P4,dT + Ke'^dba- (46) 

We write eqs.(45)-(46) in an integral form and apply the basic property of the 
Ito integral [18] 

E[ j Fdh] = 0. (47) 
Then, it follows from eq.(45) that 

3^2 

exp( — —T)E\Lr\ = const. (48) 

Eq.(48) means that the angular momentum grows exponentially in time. 

In a static metric when k = the energy po = A~^lu is also a constant of 
motion. We calculate the change of energy during the proper time diffusion (31) 
in a model without friction. Differentiating eq.(25) we obtain (see Appendix B) 

dpo = ^K^podr + KA-^Bpdh. (49) 

It follows from eqs.(49) and (47) that 

£{r) = £{s)ds (50) 

where 

£{s) = E[po{s)]. 

Eq.(50) has the solution 

f(T)=f(0)exp(^K2^) (51) 

Wc calculate a change of the angular momentum (in a spherically symmetric 
metric) resulting from the diffusion (43) with friction (Appendix B) 

^2 3^2 

dLr = —ujcpA~'^ f'LdT-l — —LrdT+KB~'^{u}-mc)p~'^Lr{pdb)+KmcB~'^xxdb. 

(52) 
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/' is negative for Jiittner as well as for quantum equilibrium distributions so that 
the angular momentum does not increase as it did in eq.(48). In the spherical 
coordinates we have 

dp^ = ^LocpA-^fdr + ^p^dr + e^dba- (53) 

ui is growing linearly with for large Pcf,. For Jiittner distribution /' = — 1 and 
for large eq.(53) is of the form 

dp^ = ap^dr - Sp'^^dr + e'^dba- (54) 

The stochastic equation for the energy is also of the form (54). We have from 
eq.(43) (see Appendix B) 

dpo = ^K^PodT + !^cl3plfdT - c(3!^^^A-'^fdT + nA-^Bpdh. (55) 

We are going to prove that the expectation value of the energy is bounded in 
time. For this purpose the Lyapunov methods [21] [22] are needed. We must 
find a non-ncgativc increasing to infinity Lyapunov function V such that for 
constants a > and K > 

gV < -KV + a. (56) 

The Ihs of eq.(56) is a time derivative of ^^[^(Pt-)]. The integral form of the 
inequality (56) reads [22] 

E[V{pr)] < exp(-/fr)F(p) + -|(1 - exp(-/fr)). (57) 
We choose V(p) — po = A~^uj. Then, by direct calculations 

Qp, = i«^Po + icPplr - c^^^A-^f. (58) 
Eq.(58) is also a consequence of eqs.(26),(47) and (55) as 

dE\pQ] = E[dpo] = QpodT 

If /' < (true for Jiittner as well as quantum equilibrium distributions) then 
from eq.(58) the inequality (56) results. As a consequence of eq.(57) the energy 
is bounded in time. The proof assumes that the stochastic equation (43) for p 
has solutions for arbitrary time. If A and are bounded functions then an 
extension of the solution to arbitrarily large time can easily be shown by means 
of the Lyapunov function method (choose V{p) = p^ as the Lyapunov function). 
We could also apply the Lyapunov argument to the stochastic equation for the 
angular momentum (53) assuming that A and are bounded ( with as 
the Lyapunov function). The problem with singular A and (the case of the 
Schwarzschild solution) is more complicated. The diffusion can terminate in the 
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singularity at finite time. The properly behaving diffusion can be constructed on 
the Kruskal extension of the Schwarzschild solution [3] . Summarizing, without 
the friction the expectation value of the energy and the angular momentum 
grow exponentially fast whereas with the friction, leading to the Jiittner or 
Bose-Einstein equilibrium distribution, the energy and the angular momentum 
are bounded in time. 

We perform some explicit (although approximate) calculations of the expec- 
tation value of the energy for the Jiittner equilibrium distribution (/' = — 1). 
Prom eq.(55) we obtain an equation 

= Jo dsSis) - 4c/3 /; dsE\po{sr] ..q. 

Eq.(59) is not a closed equation for the expectation value of the energy because 
Po is inside the expectation values on the rhs of eq.(59). We need a perturbative 
method. A linearization of eqs.(43) or (55) is not a proper tool because the Pq 
term is crucial for the equilibration. We can apply an expansion which assumes 
small K and finite k^c/S and K^vr?(? . In such an expansion the stochastic force in 
eqs.(43) and (55) can be neglected in the lowest order ( another method often 
applied in stochastic equations is based on time averaging which eliminates 
stochastic forces at the lowest order). Then, 

E{pI\ ~ E{p^f (60) 

( in general E\p^ > i?[po]^ )• Eq.(59) with the approximation (60) is equivalent 
to the differential equation 

^Sir) = (r) - ^cl3£{Tf + ^cf3K^m''c'A-^. (61) 

There is still pn in the argument of A~'^. Applying the same reasoning as 
in eq. (60) (expectation value of a function is approximately a function of the 
expectation value) we would obtain a closed integro-difFcrential equation (61) 
for £{t). Such an equation could be approached by iterative methods. We make 
a simplifying assumption that A~'^ is a slowly varying function of its argument 

£;[A-2(xO] ^A-2(x) (62) 

where the rhs does not depend on time. The solution of eq.(61) with the ap- 
proximation (62) is 

£{t) = ((£ (0) - e_)e+ eMrnHe+ - e_)) - e_ (f (0) - e+)) 

/ ^ x-i ^ (63) 

({£{0) - e_) exp(r«2(e+ - e_)) - £{0) + e+)) 

where 
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In the limit r — > cxd we obtain 



In the limit of large /3 and for the Schwarzschild solution [16] {G is the Newton 
constant, M is the mass of the source of gravity) 

^/ N 3 ._i 3 GMm 

f (00) ^ _ + mcA- ^—+mc-^. (66) 

This is the classical non-relativistic equipartition of energy in a gravitational 
potential In the limit of zero mass 

£{oo) = ^ = (^j (ipexp(-c/3|p)|) j (ip|p|exp(-c/3|p|) 

we obtain the relativistic equipartition of energy e = c£{oo) of massless particles 
(the result is correct and could be derived rigorously on the basis of an exact 
solution of eq.(80) and a calculation of its expectation values in the next section). 

The approach to the equilibrium is exponential with the speed K^(e+ — e-). 
An exponential decay to the equilibrium of relativistic diffusions coiild be proved 
by means of general methods of stochastic differential equations as discussed in 
[7], sec. 3. 2. Another method could be based on our results in [9]. We have shown 
in [9] (sec. 9) that the time evolution of the relativistic diffusion can be described 
as the time evolution in imaginary time quantum mechanics. The approach to 
the equilibrium is equivalent to the approach to the ground state (this is shown 
explicitly in the model of [10]) with the speed equal to the eigenvalue of the first 
excited state of the Hamiltonian. There are well-developed methods in quantum 
mechanics for a study of this eigenvalue. 



7 The limit m ^ 

There is a substantial simplification of stochastic equations in the limit m — >■ 0. 
We have studied this diffusion on the Minkowski space-time in [10]. We discuss 
a representation of the Poincare group for diffusing massless particles in [11]. It 
is shown that the helicity does not mix with the diffusion. Hence, our equations 
without a spin can describe diffusion of photons as well as the ultrarelativistic 
behaviour of massive particles with a spin (so it can be useful in a description of 
heavy ion collisions [23]). The limit m — >^ seems singular when applied to the 
generator (13). However, when we consider the time evolution exp(Tm^A//) 
then the limit m — >■ exists (it corresponds to a rescaling of an affine time 
parameter which in the massless case does not have the meaning of the proper 
time anyhow). The diffusion generator in the limit m — >■ reads 

Ai/=ftPfe4af^ + 3pfe^. (67) 
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The differential operator (67) defined on the manifold M does not depend on 
the metric. It is degenerate as an elliptic operator because the quadratic form 
(pjflj)^ is degenerate. As a consequence of this degeneration the diffusion in 
the configuration space is trivial. The time evolution of a massless particle 
on the Minkowski space describes a deterministic motion with a velocity of 
light on a straight line (as discussed in more detail in [20]). For this reason it 
has been rejected by Dudley [2] (sec. 11). However, in the momentum space it 
describes a diffusion of the energy. We have shown in [10] that the diffusion of 
massless particles is just a linearization of the Kompaneetz diffusion well-known 
in astrophysics [24]. On a curved manifold the motion of a massless particle 
in configuration space is non-trivial but still deterministic. The stochastic Ito 
differential equations (29)-(31) take the simple form 

— =Ai?|p|, (68) 

f = -^V (69) 

dpj = Bd^Bp^dr - cv^d^ In Adr + ^Pjdr + np^dh. (70) 

There is only one Brownian motion b for all j components of pj ( a consequence 
of the degenerate generator (67)). 
The transport equation reads 

where the friction Kk is determined by the equilibrium distribution $ e 

Kk = {i^.PjPk - IPk^E. (72) 

We discuss two examples of the equilibrium measures 

$L = Po ^ exp(-^^p0) (73) 

and 

$B =Po^exp(-/3po). (74) 

In the first case 

Kk = — ^P,pPkP<i>- (75) 
Hence, the stochastic equation for the angular momentum reads 

= -^P<t>dT - —li^pldT + Kp^db. (76) 
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Let us write 

P0 = exp(u). (77) 

The solution of eq.(76) is 



= M + K^T + Kfe^ — In ^1 + exp(u) J dsexp{K^s + Kbs)j. 
In the case (74) we obtain 



(78) 



Kk = - y/JpfcPo- (79) 
Hence, the equation for po reads 

dpo = —podr - —Ppldr + Kpodb. (80) 

It has the same solution as eq.(76) ( just replace (3^ by /3 in eq.(78)). The 
correlation functions of polynomials of po{s)~^ and p^{s)~^ can be calculated 
explicitly. The expectation values have a limit r — >■ oo expressed as moments 
of the invariant measures (73)- (74). The stochastic equations for and po do 
not depend on the metric. Their solutions and expectation values are discussed 
in [20]. 



8 A particle diffusing on de Sitter space 

In this section we discuss a relativistic particle diffusing on a background metric 
of the cosmological type 

ds^ = {dx'^f -B-^{dxf. (81) 

It is useful to apply the transformations (19)-(20) in order to reduce the stochas- 
tic equations to a perturbation of the diffusion on the Minkowski mass-shell. For 
simplicity of the formulae we restrict ourselves to B which depends only on time. 
Then, the diffusion equation (13)-(14) after the transformation (19)-(20) has the 

generator 

g = l^^m^c'Aj^^- + ±po^ - ^Bpk^ + -Lpo0B-ip,g|- (82) 
where 

Af- = (<5,fe + m-^c-^p,pk)^4^ + 3m-^c-^pk^^ (83) 

(expressed in eqs.(19)-(20) by primed momenta) is the generator of the diffusion 
on the Minkowski mass-shell. In eq. (82) 

Po = Vp^-l-m^c^ =wi. (84) 
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where uii denotes the energy on the Minkowski space-time (with the spatial 
metric tensor equal to 1). 

In the generator of the diffusion (82) the diffusion part remains the same as 
on the Minkowski mass-shell. The motion on the de Sitter space is described by 
the drift. So, the classical geodesic equations, when expressed in the coordinates 
(19)-(20), are perturbed by the relativistic Brownian motion defined on the 
mass-shell — = m^c^. 

The stochastic equations for the diffusion with the generator (82) read 

— = -uji, 85 
dr m 

— = p 86 

ctr m 

and 

dPj ^ ^cji ^^^ pjdr + ^pjdr + mcndbj + — mc)p~^pjpdh (87) 

In general, it is difficult to derive an explicit solution of such diffusion equations. 
There is one remarkable exception if (this is the metric for a causally connected 
part of de Sitter space sometimes called the Bondi-Hoyle universe [16]) 

ds'^ = idxy - cxp(— ,T°)dx2 (88) 
cR 

(R is the radius of the pscudosphere in the embedding of de Sitter space in a 
fivedimensional Minkowski space). 

Then, the drift in eq.(87) corresponding to the metric (88) is 

exactly the same as the one for the relativistic model on the Minkowski mass- 
shell (cq.(42) B = A = 1) with the friction leading to the Jiittner equilibrium 
distribution (/' = —1) and the inverse temperature 

k'/3 = ^-. (90) 

As = y^|A, we obtain a relation between the temperature and the cosmolog- 
ical constant A. Eq. (37) for the equilibrium distribution has now the solution 



$Mp)=-r^exp(-a;,^jiA) (91) 



The model of a massless particle diffusing in de Sitter space is explicitly soluble. 
Eqs.(85)-(87) read 

dx^ 
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dr " cR 

and 

dpj = --^\p\pjdT + ^-pj dr + Kpj db. (94) 

From eq.(94) it follows that Pr = jprln, where n is a time independent unit 
vector. The equation for |p| is the same as eq.(76) with the solution 

In |pt-| = In |p| + K^r + — In ^1 + -^-"2 |p| J ds exp{K^ s + nbs)^ . (95) 

We can calculate expectation values of functions of p^ using the formulae of 
Yor [25] for the probability distribution of exponentials of the Brownian motion 
and the integrals of exponentials (we have calculated some expectation values 
in [20]). The correlation functions of |Pr|~^ can be expressed by elementary 
functions. The equilibrium measure reads 

_i 2 

^ij(p) = IpI ^^p(~;^IpI)- 

It is rather surprising that the momentum evolution in de Sitter space tends to 
an equilibrium as if a certain friction was involved in a geodesic diffusion on de 
Sitter space. 



9 Summary 

We have derived basic equations concerning the relativistic diffusion with fric- 
tion in a gravitational background (adding friction to the diffusion equations of 
Franchi and Le Jan [3]) . Such equations can describe the dynamics of particles 
around stars and black holes in the presence of a gas of some other particles. 
The gravitational effects must be sufficiently strong to bo detectable by experi- 
ments. The most impressive data apply to CMBR spectrum. Then, the diffusion 
of photons scattered on charged particles in the intergalactic space can distort 
the CMBR spectrum [27]. An application of the relativistic diffusion to the 
CMBR spectrum distortion (Sunyaev-Zeldovich effect) is discussed in [10]. The 
effect of the gravitational field studied in sec. 7 is rather weak to be detected in 
a near future. Nevertheless, from the derivation of the (Kompaneetz) diffusion 
equation [24] it can be seen that if the quantum electrodynamics in a back- 
ground metric is considered then the resulting diffusion process should be the 
one perturbing the geodesic motion. The relativistic diffusion can be a useful 
tool for a selection of a relativistic approximation to multiparticle interactions 
determined by the form of the equilibrium resulting from the relativistic quan- 
tum field theory at finite temperature. We discussed as an example a diffusion 
in de Sitter space. We obtained a surprising relation between the diffusion on 
the de Sitter space and the diffusion on the flat space in a heat bath of non-zero 
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temperature. There is a well-known relation between the particle temperature 
in the de Sitter space and the radius of this space [26] . Such a relation results 
from a quantum theory and it is different from the one derived in this paper. 

The main aim of this paper is a discussion of the friction terms which lead 
to an equilibrium. It is not clear what is the physical meaning of the diffusion 
in the proper time because in contradistinction to the deterministic dynamics 
in random dynamics the proper time associated with an observer moving with 
a particle is a random variable. The transport equation (sec. 5) expressed in 
the laboratory time has a clear physical meaning. We discuss the transport 
equation in more detail in [20]. In particular, we show that the diffusion in the 
laboratory time can be obtained from the diffusion in the proper time by an 
integration over the proper time in a similar way as this is done in Feynman's 
relativistic dynamics [28] (for a discussion of the relation between the proper 
time dynamics and laboratory time dynamics see also [29]). 

10 Appendix ArStochastic equations 

In this Appendix we explain the relation between the diffusion equation (24) 

and the stochastic equations (29)-(31) (see [18] for a complete theory) . The 
main observation is that the solution at time r of the stochastic differential 
equations (assuming it is unique) is causal,i.e., depends on the initial conditions 
and values of the Brownian motion b(s) at s < t. As a consequence (of the 
Markov property) eq.(26) determines a semigroup Trcj) = 4>t- In such a case in 
order to prove that (f)T of eq.(26) is the solution of eq.(24) it is sufficient to check 
that the initial conditions are the same and the generators coincide. The initial 
condition for eq.(24) follows from the choice of the initial conditions for the 
stochastic equations. Then, the calculation of the generator is reduced to the 
solution of the stochastic equations at arbitrarily small time and the calculation 
of dcf). Let us consider eqs.(29)-(30) for and assume eq.(31) in a general form 

dpj = Bjdr + mcKejdba, (96) 

where e° are defined in eqs.(27)-(28) (together with the jjk of eq.(9) and gjk = 
B~'^5jk for isotropic coordinates). As a next step we need to calculate (for a 
small r) 

(j){xr,Vr) - (t>{x, p) ^ (t){x + Ax, p + Ap) - ^(a;, p), (97) 
where from eqs.(29)-(30) Aa;° = ^AcoAt, Ax = -^B'^pAr and from eq.(96) 

Apj = BjAr + mcKe^baiAT) - ba{0)), (98) 

here 6a (0) = 0. We expand eq.(97) in a Taylor series in Aa;^ and Apj (or what is 
the same in At and 6a (Ar)). After the expansion we calculate the expectation 
values involving the Brownian motion (see the expectation value below eq.(31); 
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we use E[ba] = and E[ba{^T)bc{AT)] = Sac^r). We must take into account 
terms till the second order in Ap which give 

E[ApjApk] = m'^c'^K^e'^elE[ba{AT)bc{AT)] = m^c^K^e^e^Ar = m^^K^^jkAT. 

(99) 

Now, collecting the terms of order At we obtain for r = At 

^A^^4>At - iB^p^.^^AT 
+5, ^ At + m-c-n- 4 afr^^- ^ ^ 

Dividing by At and taking the limit At — >^ we obtain eq.(24) with 

1 2 9^2 1 „ainA 3 2 
' 2m^ dx3 2m dx^ 2 

11 Appendix B:Time evolution of the energy 

We calculate the change of the energy 

dpo ^ d{A^^uj) = -ujA~'^djAdx^ - A^^uj^^BdjBp^dx^ + A^^B^ui^^pj o dpj. 

(101) 

Here, the circle denotes the Stratonovitch differential [18] 



hodg = hdg+ ^dhdg (102) 

The correction term to the Stratonovitch differential is 

^d{u)-'^Pj)dpj = ^m^c'^K^B-^oj-'^dT. (103) 

After an insertion of the differentials (29)-(30) and (43) on the rhs of eq.(lOl) the 
deterministic (k independent terms) cancel identically (as they should because 
Po is a constant of motion when the stochastic perturbation is absent). There 
remains 

dpo = ^A-^LJ-^B^p'^dT + Iw^^k^A-^oj-Ht + !^A-^B^cpf'p''dT 
+KA-^Bpdh. 

(104) 

This is eq.(55) (after elementary transformations). 

A proof of the formula (52) for the angular momentum is much simpler. We 

have 

rfL = rfx X p + X X dp = X X dp (105) 

We assume that the functions A and B in the isotropic metric are spherically 
symmetric. Then, in dp (defined in eq.(43)) the terms independent of k are 
parallel to x because gradients of A and B are parallel to x. The K-dependent 
terms in (105) give eq.(45). 
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